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fixed polar coordinate system. This system of equations forms the 
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ditions with inequality cons trarnH imposed on the magnitude of the 
angle of attack or the product of the dynamic pressure and the magni- 
tude of the angle of attack. The necessary conditions for optimality 
are given exclusive of derivation. Also, a computational scheme is 
given suitable for a digital computer program. 
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Symbol Definition and Units 
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reference area (m 2 ) 
vacuum specific impulse (sec) 
drag coefficient 
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launch latitude (deg) 
launch azimuth (deg) 

earth's rotational velocity in the equatorial plane 
(rad/sec) 

earth's rotational velocity in the flight plane (rad/sec) 

atmospheric pressure 

atmospheric density 

thrust (kg) 

drag force (kg) 

normal force (kg) 

gravitational acceleration (m/sec 2 ) 
gravitational constant (m 3 / sec 2 ) 
radius from center of earth to vehicle (m) 
earth-fixed velocity (m/sec) 


iv 


Symbol 


a 

« 

a 

Ro 

So 


DEFINITION OF SYMBOLS (Continued) 
Definition and Units 


earth-fixed path angle measured from local vertical 
to velocity vector (deg) 

range angle measured from earth-fixed launch point 
to radius vector (deg) 

mass of vehicle (kg-sec 2 /m)) 

angle of attack (deg) 

angle of attack rate (deg/sec) 

earth' s radius (m) 

gravitational acceleration at the surface of the 
earth (m/sec 2 ) 


M 


Mach 


X 

X* 

<p* 

™E 


attitude angle measured from the earth-fixed launch 
point to the vehicle* s thrust vector (deg) 

space-fixed attitude angle (deg) 
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earth- fixed range (m) 
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earth- fixed coordinate system 
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A SYSTEM OF EQUATIONS FOR OPTIMIZED POWERED FLIGHT TRAJECTORIES 

SUMMARY 


The equations of motion for a vehicle with thrust, lift, and drag 
forces, and a Newtonian gravitational force are derived in an earth- 
fixed polar coordinate system. This system of equations forms the dif- 
ferential equations of constraint in the calculus of variations formulation 
of minimizing flight time between two sets of boundary conditions with 
inequality constraints imposed on the magnitude of the angle of attack 
or the product of the dynamic pressure and the magnitude of the angle 
of attack. The necessary conditions for optimality are given exclusive 
of derivation. Also, a computational scheme is given suitable for a 
digital computer program. 


I. INTRODUCTION 

Normally while in the atmosphere a vehicle is constrained to fly 
a non-lifting trajectory, after tilting has been initiated shortly after 
launch, in order to minimize the structural stresses associated with 
appreciable angles of attack. In essence, this implies that the aero- 
dynamic lift is sacrificed; this may decrease the performance of the 
vehicle depending on the lift-to-drag ratio. The problem of employing 
an angle of attack during an atmospheric flight can be used if constraints 
are placed on functions related to the structural stresses on the vehicle. 
This type of problem is presented in this paper by the calculus of varia- 
tions technique where the constraint functions related to these stresses 
are explicit functions of the control variable, the angle of attack, CL. 

The purpose of this paper then is to derive the equations of motion 
in an atmospheric flight and to formulate an optimization technique with 
inequality constraints imposed on two functions related to structural 
stresses encountered during the atmospheric flight. The two functions 
chosen in this paper are the product of the dynamic pressure and the 
magnitude of the angle of attack and the magnitude of the angle of attack 
alone. 




II. DERIVATION OF THE EQUATIONS OF. MOTION 


When calculating a space vehicle’s trajectory in the atmosphere, 
it is appropriate to derive the equations of motion in a coordinate 
system fixed in (rotating with) the earth because the earth’s atmosphere 
is assumed to rotate uniformly with the earth and the external forces 
acting on the vehicle (thrust, lift and drag) are measured relative to 
the rotating atmosphere. Also, the centrifugal and Coriolis accelera- 
tion of a rotating earth affect the vehicle’s motion. 

In this paper the equations of motion are derived in a two- 
dimensional polar coordinate system. To simulate a rotating earth, it 
is assumed that this coordinate system rotates with an angular velocity 
relative to some inertial coordinate system. This angular velocity is 
measured at the launch point and is a function of the latitude and 
aiming azimuth at launch. 

In deriving the equations of motion, it is convenient to define 
two Cartesian coordinate systems. The first of these systems (denoted 
by X*) is space-fixed with its Y*-axis directed from the center of the 
earth through the launch point and with its X*-axis directed along and 
parallel to the launch azimuth. The second system (denoted by Xg) is 
earth-fixed with its coordinate axes coinciding with the space-fixed 
system at the time of launch. This system rotates with an angular 
velocity given by 


a)’ = co cos \|r sin A (2.1) 

z 

relative to the space-fixed system. In the above equation, u> is the 
earth's rotational velocity in the equatorial plane, ijr is the latitude 
above the equatorial plane, and A„ is the aiming azimuth. 

£4 

First, the space-fixed velocity vector is derived in terms of the 
earth-fixed Cartesian coordinate system. A transformation is then made 
from the earth-fixed Cartesian system to a polar coordinate system which 
is also fixed in the earth. The space-fixed acceleration vector is then 
derived in terms of the polar coordinate system. Finally, by Newton's 
second law the external forces acting on the vehicle are equated to the 
mass times the acceleration. Taking the two components of this vector 
equation yields the two scalar equations for acceleration. 
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In Figure 1, i and j are un^t vectors in the space-fixed system with 
I and j being parallel to the X and Y axis, respectively. The unit 
vectors k and 2 are parallel to the Xg and Yg axis, respectively, in 
the earth-fixed system. The unit vector e^. is defined to be parallel to 
the earth-fixed velocity vector and e^ is defined to be in the opposite 
direction of the increasing earth-fixed flight path angle, ■flg. The 
relationship between k and £ and the unit vectors in the space-fixed 
system is given by 


k = cos w' t i - sin u't j (2.2) 

~£ - sin co' t i + cos u't j. (2.3) 

Similarly, the relationship between k and £ and the unit vectors in the 
polar coordinate system is given by 

k = sin (cp E + ■flg) \ - cos (cp £ + -3 E ) e^ (2.4) 

£ = cos (cp £ + t 9 e ) e v + sin (qp £ + ■Og) e^. (2.5) 

3 


The position vector in terms of the earth-fixed Cartesian system is 


r = Xgfc + Y e I. (2.6) 

The time derivative of (2.6) in the space-fixed system gives the space- 
fixed velocity vector in terms of the earth-fixed position and velocity 
components 

? = Xglc + + Xg£ + Ygjg, (2.7) 

where the time derivatives of the unit vectors can be obtained from 
(2.2) and (2.3) which are 


k = - w* sin u't i - o)' cos u't j s ■ u' j (2.8) 

2 = co' cos u' t I - w* sin to' t 3 = u' £. (2.9) 


Substituting these values into (2.7) gives the space-fixed velocity 
vector 


r = (Xg + u' Y e ) k + (Y e - oj' X £ ) 2 . (2.10) 

The position and velocity components in (2.10) in the earth-fixed polar 
coordinate system are 


Xg = r sin cp E (2.11) 

Y = r cos <p (2.12) 

III & 

Xg = V £ sin (cp E + ^ E ) (2.13) 

Y e = V E cos (<p E + * E ). (2.14) 
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• I u 


Substituting back into (2.10) gives the space-fixed velocity vector, in 
terms of the earth-fixed polar coordinate system components, which 
becomes 


= {v E sin (cpg + -3 e ) 


+ u) f r cos 


% 


h + h 


cos (cp F + 


V - 


go r sm 9 


(2.15) 

Taking the time derivative of (2.15) in the space-fixed system and using 
(2.4) and (2.5) gives the space-fixed acceleration vector 


* S { si 


in (cp + fl ) + cos (q> + *3 ) 


^j- + id' r 


sin (cp £ + d E ) cos cp E 


- cos (<p E + -8 E ) sin cp E j- - w'r 9 £ jsin (cp E + flg) sin cp £ 


+ cos (cp E + d £ ) cos q> E J- - to' £ r -jsin (cp £ + tfg) sin 9 £ 

+ cos (cpg + -a E ) cos 9 E }] % + ["V^E + V { Sin2 (CP E + V 
+ cos 2 (9 E + d E )j- - (d'r -jsin (9 £ + -Qg) sin 9 £ 

+ cos (9 E + -9 e ) cos 9 £ j- - u'r 9 e jsin (9 £ + -a £ ) cos 9 £ 


+ cos 


(9g + ^ E ) sin C P E |' “ W ' V E |sin 2 (9 e + -Sg) + COS 2 (9 £ + -S E )j- 


- (d ,2 r -jsin (9 £ + -S £ ) cos 9 £ + cos (9 £ + -a £ ) sin 9 £ |- e^. 


(2.16) 
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In the earth-fixed polar coordinate system r and ipg are 


r = V £ cos •fig 

* V E • , 

cp_ = — sin i 9_. 
Y E r E 

Substituting these equations 
becomes 


(2.17) 

(2.18) 

into (2.16) and after simplification (2.16) 


-£ 


V— ** w* r cos 


e + 
v 


[Ve + ("T + W ' 2r ) Sin + 2 w ' Y E_ ( ’V* 


(2.19) 


This equation gives the two components of acceleration measured in the 
earth-fixed polar coordinate system. The second component is in the 
direction of the increasing flight path angle since +e^ was defined in 
the opposite direction. 



FIGURE 2 
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Assuming the vehicle as a point mass subjected to the externa] 
forces in Figure 2, the accelerations caused by these forces are equated 
to the like components in equation (2.19) by Newton's second law. 

This yields the following two scalar differential equations 


if = F cos a - - sin a - (g - w ,2 r) cos 6 

Emm i 


( 2 . 20 ) 


•8 = - " ~ sin a + 

E mV E "*’ E 


JL C0 . a + l ( t .|. rf - r ), 1 


in •8_ - 2u' , 
E 


where 


g 



(2.21) 

F - A P 

(2.22) 

v e 


I C^A - qC^A 

(2.23) 

\ P v| AO = qc^ Aa = N 4 a 

(2.24) 

GM 

7?, 

(2.25) 


where P and p are determined from some atmospheric model, ARDC or 
Patrick for example. By instantaneously considering these quantities 
as exponential functions given by 


p = P o e _ ^^ r_R °^ 


(2.26) 


p = p e^ (r_Ro) (2.27) 

K K o 9 

y and p can be obtained by taking the inverse of (2.26) and (2.27), 
respectively. The partial derivatives dP/dr and dp/ dr which are needed 
in the optimization equations would be 


dP 

dr 


7 P 


(2.28) 



(2.29) 
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Cp Q and Cjj in equations (2.23) and (2.24) are considered as functions 
of Mach number, and their derivatives are ignored in the optimization 
problem. 

The equations (2.17), (2.18), (2.19), (2.20) and the equation for 
fuel flow rate. 


F 

m = - = constant, (2.30) 

g o Is P v 

will form the system of differential equations of motion that have to 
be solved simultaneously to determine a trajectory. The next section 
will discuss the method of determining the optimum control function, 
the angle of attack as a function of time, such that a trajectory 
minimizes time between two sets of boundary conditions subject to 
equality and inequality constraints. 


III. CALCULUS OF VARIATIONS FORMULATION 


Normally while in the atmosphere a vehicle is constrained to fly 
a non-lifting trajectory, after tilting has been initiated shortly after 
launch, in order to minimize the structural stresses associated with 
appreciable angles of attack. In essence this implies that the aero- 
dynamic lift is sacrificed which may decrease the performance of the 
vehicle depending on the lift to drag ratio. The problem of employing 
an angle of attack during an atmospheric flight can be used if con- 
straints are placed on functions related to structural stresses on the 
vehicle. This type of problem can be treated by the calculus of varia- 
tions technique where the constraint functions related to these stresses 
are explicit functions of the control variable, the angle of attack, Q£. 

Two such constraint functions are used in this paper: 


g x (x, a) = QX - q|a| g 0 (3.1) 

g 2 (a) = of; - a 2 § 0, (3.2) 
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where QX and a c are numbers depending on the design limits of the 
vehicle and x denotes the set of state variables (r, Vg, -3g, cpg, m). 

The first equation (3.1) implies a constraint on the product or the 
dynamic pressure, q, and the magnitude of the angle of attack, 
whereas (3.2) constrains the magnitude of the angle of attack. 

The variational problem may now be stated as follows: It is 

desired to determine the control function, a(t), such that the expression 

tf 

M = J* f Q (x, a) dt + $>[x f , t f ] (3.3) 

t o 

is maximized subject to the inequality constraints (3.1) and (3.2) and 
the differential equations of constraint 


iL = f i (x, a, t) i = 1, ...» n 


(3.4) 


with the boundary conditions 


x. (t ) = x. 
o' 10 


(3.5) 


and 

& ( X £ j t^) i = (3.6) 

where the are the terminal constraint functions on the state 

variables. 

In order that (3.3) is maximized the following necessary conditions 
must hold (see Reference 1): 


Mk* ° 




k - 1, 2 
i ™ 1, ..., n 


(3.7) 

(3.8) 

(3.9) 
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where 


F' = f (x, a) + 


n ^ 

A i (t)|f i (x,a ) t)-x i j-+^ ^ g k (x, a), (3.10) 

frl J k=l 


and the boundary conditions 


W ■ ^ 


L v & \**t) 


(3.11) 


f o (x f> a f* + X f i^ X f ’ ^ V Z 

1=1 i-1 


(3.12) 


where the A-'s and jj^'s are Lagrangian multiplier variables and the v^'s 
are multiplier constants evaluated at t^. 

Equation (3.8) can be written as 


pF 1 _ _d_ /^F'\ 
~ dt: \^x iy / 


= 0; 


then 

N n 



(3.13) 


(3.14) 


If F' is not an explicit function of time, then (3.14) becomes 


n 



i n 


C = constant. 


(3.15) 



Then, by equations (3.8) and (3.10), equation (3.15) becomes 


n 

f o (x, a) + ^ \(t) f.(x, a) = c. (3.16) 

i=l 

If the end conditions (3.6) are not explicit functions of time, then 

c - - |£- , (3.17) 

G>t f 

since (3.16) is to be satisfied from t Q to t^. 

If we wish to minimize flight time between the boundary conditions 
(3.5) and (3.6) for a vehicle whose motion is governed by the differential 
equations of motion in part I, the expression (3.3) reduces to maximizing 


" t^ 


(3.18) 


and the function f G (x, a) is zero. Maximizing -t f is equivalent to 
minimizing t^. Equation (3.10) becomes 


F' 



(x 


z 

k=l 




(3.19) 


where the constraining equations (3.4) are given by 


V E 

m ~ sin = f x (3.20) 

f = V £ cos flg = f 2 (3.21) 

V = — — — cos a - — sin a - (g - co' 2 r) cos = f 3 (3.22) 

Emm Hi 


ii 



*E 


F - D 
mV E 


N 


V 2 

E ,2 


sin a + cos a + — ( g - — - oj' r ) sin ^ - 2w' = f 


m = - 


80 I; 


sp v 


“ f 5> 


(3.23) 

(3.24) 


and the inequality constraints are given by (3.1) and (3.2). The 
terminal constraints (3.6) of the problem are given by 


0 

II 

4 J 

s-/ 

U 

1 

u 

11 

H 

(3.25) 

* 2 - v f - v E (t f ) » 0 

(3.26) 

♦3 = - VV = °> 

(3.27) 


where rj, Vj, and #£ are the desired end conditions. 

Applying equation (3.8) to (3.19) and assuming that (3.1) and (3.2) 
are satisfied results in the following system of differential equations: 


• SF' 

A x = - = 0 

^PE 


(3.28) 


2m 


A 2 = - yjl sin a + ~ g COS ^ y (F y - F) + |jD j- cos a 


dF' _ A- 


{7 


+ mw cos <6. 


E | mV | | 7(F v " p ) + V° f sin a ~ MN cosa 
E 


! ] ' K : 


sm -3, 


(3.29) 
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^ , 2A*a 

>>=•- jiT'- A 2 cos ■ 


E L 


N sin a + D cos a 
mVi; 


(3.30) 


2Aa r . , e . 1 

j-^i j^D sin a - N cos a + sin 

[m^T ~ Sin a +H cos a + ffl|g - Y* w' 2 r^ sin J 


A4 = 


ML 


= A 2 V e sin -flg + A 3 


(u 1 r - g) sin ^ 


+ A 4 


jg(r- 8 + u ' 2r ) cos ^e] 


(3.31) 


A 5 = 


dF 

dm 


- = £(F - D) cos a - N sin aj + (F - D) sin a + N cos 


(3.32) 


Applying (3.9) to (3.14) yields the equation 


dF 1 

da 


f = - 

a 


Aa 

l v e 


A 3 (F - D + N 4 ) + A 4 §4 a si 

E j 


m a 


(3.33) 


(F - D + N 4 ) - A 3 N 4 a cos a = 0. 


The angle of attack is determined from the iterative solution of (3.33). 
If at some time the angle of attack does not satisfy the inequality 
(3.1), then a q|a| constraint exists. During this constraint period. 
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the angle of attack is computed by 


* 

a 


= qy-rg 

q a ’ 


(3.34) 


where -jr^j* determines the sign of ct. The multiplier ^ during the con- 
straint is given by 


Hi = 


L a -b=o* 


+ m 


Sg 

da 


-*0=0* 


(3.35) 


Equations (3.29) and (3.30) are augmented by 


. S£ 

' ^ ^ Sr 


Sgi Sq 

' 111 Sq * Sv E » 

respectively. The constraint ends once goes to zero. 

If the inequality (3.2) is not satisfied by the angle of attack 
from (3.33), the angle of attack becomes 


a = ± Oc, (3.36) 

where assumes the sign of a determined from (3.33). During the con- 
straint *- s given by 


M-2 = 



l- aj a=a 


-r m 


da 


J a=a 

c 


Again, the constraint ends once goes to zero. 


(3.37) 


\L 



It must be noted here that both inequality constraints cannot be 
considered simultaneously. Ifj for example, (3.1) is considered, then 
1 1 2 would be zero; or vice versa, if (3.2) is considered. 

The problem remains now of selecting the initial values of the 
Lagrangian multipliers such that all of the necessary conditions and 
the boundary conditions (3.25 - 3.27) are satisfied. First, it is 
noted that Ai(t) is a constant by equation (3.28). Furthermore, since 
cp E (t£) is allowed to vary in the problem presented in this paper, A x (t£) 
is zero by virtue of (3.11). Since f Q (x, a) is zero in equation (3.10), 
the system of equations (3.29 - 3.32) is homogeneous in the A's. This 
property makes one of the initial values of the A's arbitrary. In this 
paper, we have chosen 


A3 (to) _ A30 58 1 • 


(3.38) 


By specifying an initial angle of attack, 0 ^, and using (3.33), the 
initial value of A* is given by 



(3.39) 


Furthermore, by specifying an initial angle of attack rate, Qq, and 
differentiating (3.33) with respect to time and using (3.28), (3.29), 

A 3 o> and A 4 o> A 20 can be determined. Since equations (3.29 - 3.32) 
are independent of A 5 , the initial value of A 5 can be taken as zero. 

Thus, Oo a nd cto a re used to isolate the terminal constraints (3.25 - 
3.27). Once these terminal constraints have been satisfied the terminal 
values of the A's are by (3.11) and (3.25 - 3.27) 


Ag(tf ) = Vi, A3(tp - V2> and A^t^) - V3. 

Since, in equation (3.12) the ijr/s are not explicit functions of time, 
(3.12) is given by 


n 

Y A ± (t f ) f£(x f , a f )= - = const, (3.40) 

i=l f 

which becomes (3.16) evaluated at t£. 
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IV. APPLICATIONS FOR A DIGITAL COMPUTER PROGRAM 


In this section, a computational scheme is given suitable for a 
digital computer program. First, the inequality constraints on q 1 0i| 
and a are ignored. This is the basic program, and only slight modifica- 
tions are needed to include the inequality constraints which are given 
after the basic scheme. 

A. Application Without Inequality Constraints 
1. Input Data Needed 

Constants: g , R , C x , C 2 , A , i|/, to, At, 

o o z 


Initial Conditions: cp„ , r , V , •a 

t 0 o t 0 b 0 


Tables: C Dq , C^, Mach 


Isolation Parameters: OLq, OCq 


Vehicle Data: F^Clbs), 

I S p v (sec), A (m 2 ) , A g (m 2 ), W Q (lbs) 

Pr e 1 oad Compu ta t ions 


o 9 

II 

O 

+ 

n 

H 

(4.1) 

F 


. v 

m ” C I 

1 Sp y 

(4.2) 

F v (kg) = F v (lbs) -s- C 2 

(4.3) 

co* = oj sin A z cos y 

(4.4) 

GM = g 0 R§. 

(4.5) 
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To compute A 20 an ^ ^ 40 * l et 


doc 


„ . 

K * ' " Ha ’ 


Then equation (3.33) becomes 


f a = A 3 j 4 + A 4 K 4 = 0 


at t Q : A 30 = 1, a = Ob- 


Then, 


A40 “ ” R 4 4- J4 


Taking the time derivative of (4.8) gives 


A 3 J 4 + A 3 J 4 + A 4 K 4 + A 4 K 4 * 0. 


To determine J 4 and K 4 , we compute: 


„= — In ~ 

^ r ” R o P 

1 , fo 

7 " r - R„ n P 


Dl . a - . ^ 


2 av E v E 


(4.6) 

(4.7) 

(4.8) 


(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 
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(4.15) 


N 



|oN 


„ _ SN _ 2N 
2 = ^ V E " V E 


M _ M± - KF 

N 41 “ ” " M^4 


N 


42 


_ ^ 
av r 


2N, 

V„ 


N = N 4 a + N^r + N 2 V e 
D - Dj.* + D 2 V e 
N 4 = N 41 f + N 42 V e 


K 4 i = 





cos a + Ni 


sm a 


(4.16) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 


K 


_ = _1_ 

42 * SV E ‘ mV E 


(D - N 42 ) cos a + N 2 sin a - mK 4 


(4.23) 


K. 


= ^ m _1_ 
44 da ~ mV_ 


|(F - D) + 2N 4 | sin a + N cos a 


(4.24) 


then 




- mJ 4 + *|n 4 + yx (F y - F) 


- d - Naj- 


sm a 


+ -i (N 4 + F - D) a + N y cos a 


(4.25) 


1 ft 



and 


K 4 = K44Q! + K 41 r + K 42 Vg - ~ K 4 . ( 4 . 26 ) 

Substituting equations (3.29) and (3.30) for A 3 and A 4 into (4.10) and 
evaluating at t D with A 30 and A 40 and using the above equations A 2 o 
becomes 


A 2 o ~ 4 - E 2 , 


(4.27) 


where 


L mV E, 


(N sin Oo + D cos Ofo) 


” "^4 d" A30 E4 + J4 

mV 2 

*£ ^ 

+ ( D 8in Oo - N cos % + Sin ^ J + 

+ &-}]+ K 4 [(m'* r G - 8) sin ^ 




_ s _ + £lL 
X X 


') “ 8 xj 


(4.28) 


and 


E 2 * J 4 cos K 4 Vg sin ^g 1 

a o ° o 


(4.29) 


3. Trajectory Integration and Isolation 

At a time, t N , the equations of motion (3.20 - 3.24) and 
the "lambda dot" equations (3.20 - 3.32) are numerically integrated 
from t N to ttffi, where the variables r, V £ , -dg, A 2 , X 3 , A 4 , and A 5 
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needed to start the integration are given, have been precomputed, or 
have been integrated from the t N-1 time point. The atmospheric pres- 
sure and density are determined from an atmospheric model subroutine 
incorporated into the program. The lift and drag coefficients, and 
C Do , are determined by Lagrangian interpolation as functions of Mach 
number. The angle of attack at t N is determined from the following 
subroutine. 


where 


Compute f a : The initial a to begin the iteration is taken 

from the fc N- time point. 

If f^, g Tol., a has been determined. (Tol = Tolerance) 

If f a > Tol., we go to equations (4.30): 



(4.30) 


df 


f 1 = 

a 


ha 


a _ -v 

- a 3 


(F - D + 2N 4 ) cos a - N sin a 



D + 2N 4 ) s in a + N cos a . 


(4.31) 


f 

If g Tol, then a has been determined. 

a 


If j-jrj > Tol, compute a new a from (4.30) using the 

a 

preceding a to compute f^, and f^, , and repeat this pro- 
f 

cedure until j-jr| g Tol. 

a 
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Trajectory Isolation; 


A trajectory cuts off on local space-fixed circular 
velocity - the program may be modified to cut off on any velocity - and 
isolates altitude and the space-fixed flight path angle. Since the 
vehicle's motion is measured in an earth-fixed coordinate system, a 
transformation is used to compute the space-fixed values for velocity 
and the path angle. The relationships are 


V* = V 2 + 2i o' 


-, 1 / 2 


Vj,r sin •Sg + w^ 2 cos 2 \|r 


and 


43* 


tan -1 


Vg/V* cos 43g 
Wl - (Vg/V*) 2 cos 2 flg' - 


(4.33) 


(4.34) 


(Xq and do are then used to isolate these end conditions. The isolation 
scheme given in Reference 2 is well suited for this problem. 

4. Additional Equations Useful in Trajectory Analysis 

The attitude of the vehicle measured from the vertical 
earth-fixed launch point is given by 

x = cp E + ^ + a. (4.35) 

and the space-fixed attitude angle is given by 

X* = X + U)'t. ^ (4.36) 


The range of the vehicle measured on the surface of the earth from the 
earth-fixed launch point is given by 


XXX = R 0 cpg. 


(4.37) 
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and the range measured from the space-fixed launch point is 


XXX* = XXX + R 0 U)' t . 


(4.38) 


B. Application with qlq] Constraint 

Input Data (In Addition to IIIA) : QX, Ago* A 40 Here QX is 

a constant which q|a| cannot exceed during the trajectory; A 20 and A 40 
are used as the isolation parameters rather than Oo and Oo, since the 
angle of attack is determined from the constraint equation and cannot 
be used as an isolation parameter. The computation of equations (4.9 - 
4.29) is ignored in this modification of the program. 

The following subroutine is used: 

We compute: 

81 = QX - q|a|, (4.39) 

where a has been determined from the iteration of f^. 

If gx § 0, Hi = 0 and the angle of attack determined from the 
f equation are used in the main routine of the program. If g x < 0, 
then the angle of attack becomes 


and 




a 


q a 


m - 




+ m 



(4.40) 


(4.41) 


Equations (3.29) and (3.30) are then augmented by 


- hi (nq|«|) 

hi (pv e M>., 

respectively. The constraint period ends when hi goes to zero. 
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C. Angle of Attack Constraint 
Input needed ; o^, A 2 o> A40 
The following subroutine is used: 

We compute: 

g 2 = " a 2 . (4.42) 

If g £ g 0, Hi = 0 and the angle of attack determined from the 
f a equation are used in the main routine of the program. 

If g 2 < 0, then the angle of attack becomes 


a = ± etc* 


(4.43) 


where assumes the sign of a determined from the f a equation and 


H2 " 




(4.44) 


The constraint period ends where |i 2 goes to zero. 


V. CONCLUSIONS 

The optimization technique given in this report has been success- 
fully applied to the Saturn V and Saturn IB vehicles. The equations 
were programmed for the IBM 7094 digital computer. The results, which 
are to be published later, show that an increase in orbital payload 
can be obtained as compared to the non-lifting first stage trajectories. 
This increase of orbital payload is consistent when the q|o:| product is 
constrained to a reasonable value. The existing program does not con- 
tain the angle of attack constraint at the present time, but it is felt 
that this can be easily incorporated into the program. 
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